A generic new platform for topological quantum computation using semiconductor heterostructures 
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We show that a film of a semiconductor in which s-wave superconductivity and a Zeeman sphtting are in- 
duced by proximity effect, supports zero-energy Majorana fermion modes in the ordinary vortex excitations. 
Since time reversal symmetry is explicitly broken, the edge of the film constitutes a chiral Majorana wire. The 
heterostructure we propose - a semiconducting thin film sandwiched between an s-wave superconductor and a 
magnetic insulator - is a generic system which can be used as the platform for topological quantum computation 
by virtue of the existence of non-Abelian Majorana fermions. 
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Introduction. In two spatial dimensions, where permuta- 
tion and exchange are not necessarily equivalent, particles can 
have quantum statistics which are strikingly different from 
the familiar statistics of bosons and fermions. In situations 
where the many body ground state wave-function is a lin- 
ear combination of states from a degenerate subspace, a pair- 
wise exchange of the particle coordinates can unitarily ro- 
tate the ground state wave-function in the degenerate sub- 
space. In this case, the exchange statistics is given by a multi- 
dimensional unitary matrix representation (as opposed to just 
a phase factor) of the 2D braid group, and the statistics is non- 
Abelian [1]. It has been proposed that such a system, where 
the ground state degeneracy is protected by a gap from lo- 
cal perturbations, can be used as a fault- tolerant platform for 
topological quantum computation (TQC) 1 2] . 

Recently, the u = 5/2 FQH state at high magnetic fields 
and at low temperature has been proposed as a topological 
qubit |0t]. This theoretical conjecture, however, awaits ex- 
perimental verification ||3t |4t] . An equivalent system, in which 
the ordered state is in the same universality class as the 5/2 
FQH state, is the spin-less (spin-polarized) Px + ipy super- 
conductor/superfluid |5t]. In a finite magnetic field, a vor- 
tex excitation in such a superconductor traps a single, non- 
degenerate, zero-energy bound state. The key to non-Abelian 
statistics is that the second-quantized operator for this zero- 
energy state is self-hermitian, 7^ = 7, rendering 7 a Majorana 
fermion operator. If the constituent fermions have spin, the 
spin-degeneracy of the zero energy excitation spoils the non- 
Abelian statistics, because the mutual statistics of the vortex- 
fermion composites becomes trivial. To circumvent this prob- 
lem in a realistic, spinful, Px + ipy superconductor such as 
strontium ruthenate, it has been proposed that the requisite 
excitations are the exotic half-quantum vortices, which can be 
thought of as ordinary vortex excitations in only one of the 
spin sector in the condensate [6] . 

Even though quenching the spin-degeneracy by either the 
application of a magnetic field ^ or by using spin-less 
atomic systems fs] is possible in principle, it is practically 
very difficult. Therefore, it is desirable to have systems whose 
most natural excitations themselves follow non-Abelian statis- 
tics in spite of the fermions carrying, as they do in a realistic 



system, a spin quantum number. The recent proposal by Fu 
and Kane [9] points out one such system - the surface of a 
strong TI in proximity to an 5-wave superconductor - which 
supports a non-degenerate Majorana fermion excitation in the 
core of an ordinary vortex. In this paper, we propose a simple 
generic TQC platform by showing that it is possible to re- 
place the TI with a regular semiconductor film with spin-orbit 
coupling, provided the time-reversal symmetry is broken by 
proximity of the film to a magnetic insulator. The three ingre- 
dients of non-Abelian statistics - spin-orbit coupling, s-wave 
superconductivity, and Zeeman splitting - are experimentally 
known to occur in many solid state materials. It is encourag- 
ing that the 5-wave proximity effect has already been demon- 
strated in 2D InAs heterostructures which additionally also 
have a substantial spin-orbit coupling Thus, the struc- 
ture we propose is one of the simplest to realize non-Abelian 
Majorana fermions in the solid-state context. 

Theoretical Model. The single-particle effective Hamilto- 
nian for the conduction band of a spin-orbit coupled semi- 
conductor in contact with a magnetic insulator is given by (we 
set = 1 henceforth) 

HQ = —--ii^VzCrz^a{axp)-z. (1) 
2m* 

Here, m*, Vz and /i are the conduction-band effective mass of 
an electron, effective Zeeman coupling induced by proximity 
to a magnetic insulator (we neglect the direct coupling of the 
electrons with the magnetic field from the magnetic insulator), 
and chemical potential, respectively. The coefficient a de- 
scribes the strength of the Rashba spin-orbit coupling and cr^ 
are the Pauli matrices. Despite the similarity in the spin-orbit- 
coupling terms, Hq and the Hamiltonian for the TI surface in 
Ref . [£] differ by the existence of a spin-diagonal kinetic en- 
ergy term in Eq. O. Because of the spin-diagonal kinetic en- 
ergy, there are in general two spin-orbit- split Fermi surfaces in 
the present system, in contrast to the surface of a TI in which 
an odd number of bands cross the Fermi level ii. In Eq. O, 
for out-of-plane Zeeman coupling such that \ Vz\ > a sin- 
gle band crosses the Fermi level. Thus, analogous to a strong 
TI surface (but arising from qualitatively different physics), 
the system has a single Fermi surface, which is suggestive of 



2 




Magnetic insulator 



FIG. 1 : (Color online) Schematic picture of the proposed heterostruc- 
ture exhibiting Majorana zero-energy bound state inside an ordinary 
vortex. 



non-Abelian topological order if s-wave superconductivity is 
induced in the film. We show below that this is indeed the case 
by analyzing the Bogoliubov de Gennes (BdG) equations for 
a vortex in the superconductor in the hetero structure shown in 

Fig.m 

The proximity-induced superconductivity in the semicon- 
ductor can be described by the Hamiltonian, 



I dr{Ao(r)c|(r)c|(r)+H.c}, (2) 



where cj. (r) are the creation operators for electrons with spin 
a and Ao(r) is the proximity-induced gap. The corresponding 
BdG equations written in Nambu space become, 
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where ^(r) is the wave function in the Nambu spinor basis, 
^(r) = (ii^(r), 'U|(r), 'U|(r), — 'u^(r))^. Using the solutions 
of the BdG equations, one can define Bogoliubov quasiparti- 
cle operators as 7 = J dr '^^Ua{r)c'l{r)-\-Va{T)ca{T). The 
bulk excitation spectrum of the BdG equations with A(r) = 
Ao has a gap for non- vanishing spin-orbit coupling. 

BdG equations for a vortex. We now consider the vortex 
in the hetero structure shown in FigH] and take the vortex-like 
configuration of the order parameter: Aq (r,6>) = Ao(r)e*^. 
Because of the rotational symmetry, the BdG equations can 
be decoupled into angular momentum channels indexed by m 
with the corresponding spinor wave-function, 

*„M) = e™^ {u^{r),u^{r)e'',v^{r)e-'',-v^{r)f. (4) 

Because of the particle-hole symmetry of the BdG equations, 
if ^(r) is a solution with energy E then iayry'^'^{r) is also 
a solution at energy —E. Here Ty is defined to be the Pauli 
matrix in Nambu spinor space. In particular, zero-energy 
solutions of the BdG equations must come in pairs, ^(r) 
and iayry'^'^{r), unless these two wave-functions refer to 
the same state. Thus, a zero-energy solution in an angular 
momentum channel m is always paired with another zero- 
energy solution in the channel —m and, therefore, can be non- 
degenerate only if it corresponds to the m = angular mo- 
mentum channel. 



The radial BdG equations describing the zero-energy state 
in the m = channel can be written as 
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^(r) = 0, 



(5) 
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with T] = Since the BdG matrix is real, there are two so- 
lutions ^(r) and ^*(r) with E = 0. Furthermore, it follows 
from the particle-hole symmetry of the BdG equations that 
i(jyTy^{r) is also a solution. Thus, any non-degenerate E = 
solution must satisfy the property i(jyTy^{r) = zA^(r). 
Moreover, because {icFyTy)'^ = — 1, the possible values of A 
are A = ±1. The value of A sets a constraint on the spin- 
degree of freedom of ^(r), such that v^{r) = \u^{r) and 
u^{r) — \v^(t). This allows one to eliminate the spin de- 
grees of freedom in ^ (r) and define a reduced spinor (^) = 
(ix^(r), ?x^(r))^. The corresponding reduced BdG equations 
take the form of a 2 x 2 matrix differential equation: 



(-r^(dl^\dMV,-^ AA(r)+a(S, + ^) 



^o(r) = 0. 



(6) 



We now approximate the radial dependence of Ao(r) by 
Ao(r) = for r < and Ao(r) = Aq for r > i?. In 
view of the stability of the putative Majorana zero-energy so- 
lution to local changes in the Hamiltonian [5], such an ap- 
proximation can be made without loss of generality. For 
r < the analytical solution of Eq. © is given by ^o(^) = 
(u^Jo^zr\ uiJi{zr))^ with the constraint 
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0. (7) 



Here Jn(r) are the Bessel functions of the first kind. The 
characteristic equation for z is 



ivz'-^y 



a^z'^ = 0. 



(8) 



In the case /j. > Vz the roots of Eq.® are real: zi^2 = 
^-s/wi and ^3,4 = with wi^2 > being the solu- 

tions for z^. In the opposite limit, < /j. < Vz, there 
are two real solutions zi^2 — ■^\l\w\\ and two imaginary 
solutions ^3,4 = ±.^^J\w2\• Since the Bessel functions are 
symmetric, we find two solutions which are well-behaved at 
the origin: ^i(r) = {u^Jo[z\t),u^J\{z\t)Y and ^2(^) = 
{u^Jo{^zzr\u^J\{^zzryj^ . Therefore, the full solution at r < 
i?is \E^o (0 = ci^i(r) +C2^2(r). 

At large distances r > R, where Ao(r) = Aq, the solu- 
tion to Eq. S is complicated. Nevertheless, one can write the 
solution as a series expansion in 1/r: 
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FIG. 2: (Color online) Upper panel: complex roots of the Eq. ([Tol) for 
different values of /i and A with Ao / 0. Lower panel: Numerical 
solution for the Majorana zero-energy state ^^{r) — {u{r)^v{r))^ 
for A = —1 and fi < Vz. The dashed (red) and solid (blue) lines 
correspond to u(r) and v{r), respectively. Here we used the fol- 
lowing parameters: rj — a — Vz — l,/i = 0, Aq = 0.1 and 
R = 1. The boundary conditions used are ^o(O) = (IjO)^ and 
^o(r = 40) = (0,0)^. 

where are the corresponding spinors. The zeroth order co- 
efficient ao satisfies the following equation: 

y — AAo — iza r]z^ — Vz — jJ^ J 

The higher order coefficients can be calculated from ao 
using a set of recursion relations |11]. The characteristic 
equation for Eq. ([TOb has 4 complex roots for z, which are 
shown in Fig [21 Physical solutions of Eq. ([6l) at r > R, 
^o(^) = Z]n>2^n^n(r), require that Im[2;n] > 0. (Here 
^n(^) is the solution corresponding to the eigenvalue Zn) 
Thus, for (/i^ + ^o) ^ K^' ther^ are two solutions for 
A = ±1. On the other hand, for (/i^ + Ag) < there 
are three solutions for A = —1 and only one for A = 1. 

In order to obtain a unique solution for the zero-energy 
state, the wavefunctions (^) (^) should satisfy 

boundary conditions aX r = R. Since we are matching 2- 
component wave-functions and their derivatives air = R, the 
continuity of ^o(i?) and ^o(i^) leads to 4 independent equa- 
tions. One additional constraint comes from the normalization 
of the wavef unction in all space. Thus, there are five indepen- 
dent constraints for the coefficients c^. A unique zero-energy 
solution exists if the number of unknown coefficients Cn is 
five, which is the case for (/i^ + Aq) < and A = —1. In 
this case, the wavefunctions (^) — ^n=i 2 Cn^n(^) and 
^0 (^) — Sn=3 4 5 <^n^n(^) havc 2 and 3 unknown coeffi- 
cients, respectively. In all other cases the number of unknowns 
Cn is smaller than five, and thus, as we have checked explicitly, 
solutions for the zero-energy eigenfunction do not exist. From 
these arguments one can conclude that, for (/i^ + Aq) < V^, 
an ordinary vortex in the superconducting condensate contains 
a unique non-degenerate E = solution in the m = angu- 



lar momentum channel. The numerical solution for the zero- 
energy state is shown in Fig. [2^. It is straightforward to check 
that the zero-energy solution corresponds to a self-hermitian 
second-quantized operator 7 = 7^: it is a Majorana fermion 
excitation. 

We can also consider the special case with r] = O^Vz = in 
the above equations, which describes the recent proposal for 
TQC [9] using zero-energy Majorana bound states at vortices 
on the interface of a TI and an s-wave superconductor. In this 
case, we find a single solution fovr < R and a pair of indepen- 
dent solutions for r > R. Since the BdG differential equation 
is now only first order, we need only match the 2-component 
spinors themselves (derivatives need not match) which yields 
3 equations for the 3 coefficients. This leads to a unique 
Majarana fermion solution at the vortex, which is consistent 
with Ref . 1 9] . Interestingly, in contrast to our Hamiltonian for 
?7 > 0, the condition for the existence of a Majorana fermion 
for 7^ = is given by < (Aq + /i^). The model con- 
sidered in Ref. [9] and our present system have a similar order 
parameter structure. In both cases the order parameter compo- 
nent {c^{r)ci{r')) has an 5- wave orbital symmetry while the 
order parameter components {c^{r)c^{r')) and {ci{r)ci{r')) 
have px + ipy and px — ipy orbital symmetries, respectively. 
On the surface of a TI, because of time-reversal invariance, 
\{c^{r)c^{r^))\ = |(c|(r)c|(r'))|. In our system, the ratio of 
the order parameter components in the two spin sectors is dif- 
ferent from 1, and approaches 1 in the limit a'^/r] Vz. In 
both cases, however, the superconducting pairing potential is 
s-wave, and is induced by proximity effect. Therefore, the 
superconducting state and the associated non-Abelian topo- 
logical character are expected to be robust in the presence of 
finite disorder. 

Topological phase transition. We have shown above that 
a non-degenerate Majorana state exists in a vortex in the su- 
perconductor only in the parameter regime (/i^ + Aq) < V^. 
This suggests that there must be a quantum phase transition 
(QPT) separating the parameter regimes (/i^ + Aq) < V^ and 
(/i^ + Aq ) > V^, even though the system in both regimes is an 
5-wave superconductor. A non-degenerate zero-energy solu- 
tion cannot disappear unless a continuum of energy levels ap- 
pears around ^ = 0. Such a continuum of states at £^ = can 
only appear if the bulk gap closes, which can be used to define 
a topological quantum phase transition. In the present system, 
such a phase transition can be accessed by varying either the 
Zeeman splitting or the chemical potential. A similar topolog- 
ical quantum phase transition has already been predicted for 
ultra-cold atoms with vortices in the spin-orbit coupling 1121 . 

The bulk gap of the present system can be calculated from 
the bulk excitation spectrum, 

= V^ + Ag + g2 ^ a^k^ ± 2^V^^l^l\V^^OL^k^) 

(11) 

where e = r]k'^ — /j.. As seen in Fig. [3l the excitation gap 
first increases as a function of Aq (proximity-induced pair- 
potential) and then decreases and vanishes at a critical point, 
Aqc = \/V^ — jj?, before re-opening and increasing with 
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FIG. 3: (Color online) Quasiparticle gap versus pairing potential for 
various values of the chemical potential /i. Here Eg^^ed = Eg/Vz, 
Ao,red = A0/V2, /ired = ^7^2 and ttred = a/V^T^. Thc max- 
imum value of the gap in the topologically non-trivial superconduc- 
tor, and the corresponding area in the phase diagram, decreases with 
increasing values of /x. For large negative /i, the system makes a 
transition to a semiconductor. The phase to the right of the critical 
point is the topologically trivial s-wave superconductor. 



Aq. The critical point marks the phase transition between 
a topologically non-trivial (left) and a topologically trivial 
(right) s-wave superconducting phases. The scale of the gap 
in the topologically non-trivial phase is set by the strength of 
the spin-orbit coupling a and the position of the critical point. 
The fact that the phase on the right-side of the critical point 
does not support a non-degenerate Majorana mode can be ver- 
ified by observing that, for these values of Aq, it is possible 
to reduce Vz such that \ Vz\ < without the gap vanishing at 
any point. This is the phase without Majorana Fermion excita- 
tions. In fact, this phase can be reached from the conventional 
5-wave superconductor with Vz = and a = without cross- 
ing a phase transition. 

Majorana edge modes and TQC. In analogy with Ref . , 
we find that an interface between two superconductor layers, 
which can be deposited on the semiconductor thin film, sup- 
ports a pair of zero-energy excitations when the phase differ- 
ence between the superconductors is tt. This geometry can be 
analyzed in a way that closely follows our derivation of the lo- 
calized state in a vortex in the m = channel, since the BdG 
Hamiltonian can again be reduced to a real matrix. In this 
case, we find that, in the parameter regime (/i^ + Aq) < V^, 
there are 3 linearly independent solutions on each side of 
the interface. Since the number of constraints to be satisfied 
at the interface (we assume the interface to be of negligible 
width) remains 5 as before, one expects a pair of indepen- 
dent zero-energy solutions. The interface, therefore, consti- 
tutes a non-chiral Majorana wire, which can be exploited for 
braiding in a way completely analogous to Ref. to perform 
TQC. Majorana bound states as well as Majorana edge modes 
in our system can be studied experimentally using non-local 
Andreev reflection [13] and electrically detected interferome- 
try UAilfJ experiments. 



The experimental implementation of this proposal involves 
a heterostructure of a magnetic insulator {e.g. EuO), a strong 
spin-orbit coupled semiconductor {e.g. In As) and an s-wave 
superconductor with a large Tc {e.g. Nb). Using these mate- 
rials, it is possible [H] to induce an effective superconduct- 
ing pairing potential Aq ~ 0.5 meV and a tunneling-induced 
effective Zeeman splitting Vz ^ I meV. Additionally, the 
strength of spin-orbit interaction a in InAs heterostructures 
is electric-field tunable and can be made as large as a ^ 50 
meV-A [16]. With these estimates, the quasiparticle gap Eg is 
of the order of 1 K. Given that the chemical potential is gate- 
tunable and can be of the of the order of Aq, we numerically 
estimate the magnitude of the excitation energy for the bound 
states in a vortex core of size ~ 20 nm to be of the order of 



0.1 K lUlh . which sets the temperature scale for TQC in this 
system. 

Conclusion. Our proposed TQC platform should be sim- 
pler to implement experimentally than any of the TQC candi- 
dates proposed in the literature so far, since it involves a stan- 
dard heterostructure with a magnetic insulator, a semiconduc- 
tor film, and an ordinary 5-wave superconductor. We believe 
that the proposed scheme provides the most straightforward 
method for the solid-state realization of non-Abelian Majo- 
rana fermions. A significant practical advantage of the pro- 
posed TQC scheme is its generic simplicity: it requires nei- 
ther special samples or materials nor ultra-low temperatures 
or high magnetic fields. 
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